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I.     ABSTRACT 

Strong  rain-induced  mixing  in  a  thin  surface  layer  is  numerically  shown  to  greatly 
increase  surface  gravity  wave  attenuation.  This  case  study  uses  a  single  wavelength  (2.8 
m)  together  with  two  mixed  layer  depths  (10  and  20  cm).  The  rain-induced  mixing  is 
simulated  by  varying  kinematic  viscosity  within  the  mixed  layer  from  10"*  to  10"^  m-  r*'  , 
molecular  to  strong  turbulent  mixing,  respectively.  The  results  indicate  that  surface 
gravity  wave  attenuation  in  the  presence  of  a  thin  rain-induced  mixed  layer  can  increase 
by  a  factor  of  up  to  60U0  times  the  attenuation  rate  due  to  molecular  viscosity  alone. 
This  indicates  that  rain  need  only  mix  the  top  10-20  cm  surface  layer  to  effectively 
dampen  short  surface  gravity  waves. 


II.     INTRODUCTION 

For  centuries  sailors  and  casual  observers  living  along  the  coast  have  observed  that 
ocean  waves  are  dampened,  and  the  sea  becomes  flat  during  rain.  This  dampening  is 
evident  especially  during  rain  events  when  winds  are  light  or  calm.  This  dampening  must 
be  a  result  of  energy  transfer  from  a  raindrop  that  strikes  the  ocean  surface.  The 
raindrop  impact  may  generate  capillary  waves,  subsurface  turbulent  mixing,  and  sound 
or  heat.  Only  subsurface  turbulent  mixing  is  likely  to  interact  with  surface  gravity 
waves.  This  process,  first  suggested  by  Re}nolds  (1900),  takes  the  form  of  vortex  rings 
which  propagate  downward,  mixing  the  water  in  a  near  surface  layer  as  the  rings  disin- 
tegrate; thus  enabling  the  vortex  rings  to  eflectively  increase  the  turbulent  mixing  in  a 
thin  layer  near  the  surface  of  the  water. 

Manton  (19' 3)  used  an  analytical  approach  based  on  Reynolds'  idea,  to  argue  that 
rain-induced  mixing  will  cause  surface  gravity  wave  attenuation.  Manton  assumed  a 
constant  eddy  viscosity  throughout  the  water  column  to  solve  his  governing  momentum 
equations.  The  assumption  of  no  depth  dependence  of  the  mixing  seems  unreahstic. 
The  scope  of  this  study  is  to  examine,  by  numerical  methods,  the  effects  of  increased 
turbulent  mixing  due  to  rain  in  a  thin  surface  layer.  It  is  hypothesized  that  by  increasing 
the  turbulent  mixing  in  a  thin  surface  layer,  significant  attenuation  of  a  surface  gravity 
wave  will  result.  This  study  focuses  on  waves  having  a  wavelength  of  2.8  m  and  ampli- 
tudes of  3  and  6  cm.  Eddy  viscosity  is  varied  from  10"^  to  10"^  tri^  5~'  .  The  mixed  layer 
depth  is  either  10  cm  or  20  cm.  The  rate  of  kinetic  energy  change  at  each  time  step  was 
calculated  and  used  to  determine  the  rate  of  wave  attenuation.  The  use  of  a  numerical 
calculation  allowed  viscosity  to  be  varied  with  depth.    Altcrnativch .  Manton  employed 


a  velocity  potential  solution;  however,  this  precluded  solving  for  the  rotational  flous 
occurring  here. 

A.     SIGNIFICANCE 

The  dampening  of  waves  by  rain  can  influence  sateUite  remote  sensing  observations 
at  the  sea  surface,  in  addition  to  changing  the  characteristics  of  the  water  column.  Sat- 
ellite borne  scatterometers  and  Synthetic-Aperture  Radars  (SAR)  measure  microwave 
energy  backscattered  to  space  from  the  ocean's  surface.  Scatterometers  are  able  to  de- 
termine wind  speeds  based  on  the  energy  backscattered  from  the  sea  surface  by  short 
surface  gravity  waves  (approximately  30  cm).  If  these  waves  are  attenuated  by  rain, 
scatterometer-derived  winds  will  be  in  error  for  areas  under  a  rain  event,  giving  the  im- 
pression that  there  are  light  winds  present  when  in  fact  greater  winds  may  be  present. 

Fu  and  Holt  (1^82)  noticed  that  squall  storm  lines  were  easily  detected  in  SAR  im- 
ages. They  found  that  down  draft  and  rain  areas  were  delineated  by  changes  in  image 
brightness.  The  image  was  brighter  in  the  outflow  area  due  to  a  rough  sea  surface.  A 
smooth  dark  image  caused  by  tlie  absence  of  backscattered  energy  suggested  that  m 
areas  of  rain  the  capillary  and  short  gravity  waves  had  been  dampened. 

Rain  also  iniluences  tiie  water  column.  Livezey  (19SS)  found  that  rain  events  could 
cause  changes  to  the  sea^^onal  mixed  layer  depth  indirectly  effecting  heat  loss  to  the  at- 
mosphere. Using  a  total  kinetic  energy  budget  model.  Garwood  (1989)  demonstrated 
that  a  rain-generated  mixed  layer  would  form  in  a  stratified  water  column.  .-\  2'>()  cm 
mixed  layer  depth  (MLD)  was  produced  when  turbulent  kinetic  energy  (TKE)  of  a  rain 
event  was  injected  into  the  upper  layer  of  the  ocean,  see  Fig.  1.  Within  1(X»0  seconds 
the  MLD  is  seen  to  be  approaching  steady  state.  A  continued  deepening  of  the  MLD 
will  occur,  but  at  a  much  slower  rate  until  equilibrium  has  been  reached.  The  results  are 
tentative,  but  do  suggest  that  rain  has  a  possible  role  in  ocean  mixed  layer  dynamics. 
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Figure   I.       Mixed  layer  depth  prediction:     A  200  cm  mixed  layer  has  formed  by 
rain-injected  TKE. 


When  a  raindrop  strikes  a  fluid  surface  there  are  other  interesting  efTects.  One  effect 
is  a  roughening  of  the  sea  surface  from  depression  cusps  that  break  the  smoothness  of 
the  sea  surface.  In  conjunction  with  the  wind,  the  roughened  surface  caused  an  increase 
in  the  surface  stress.  This  stress  increase  augments  the  transfer  of  energy  into  the  water 
column.  Caldwell  and  Elliott  (1970)  found  that  a  rain  rate  of  several  centimeters  per 
hour  in  light  winds  could  produce  stresses  comparable  in  magnitude  to  the  wind  stress 


for  winds  as  high  as  10-20  m  ^'  without  rain.  A  heavy  rain  can  also  produce  bubble 
clouds,  a  layer  of  bubbles  beneath  the  wave  surface,  similar  to  those  produced  by 
breaking  waves.  Thorpe  (1984)  used  these  bubble  clouds  to  examine  the  effects  of 
Langmuir  circulation  by  following  their  progress  through  the  water  column.  In  his  study 
he  found  that  the  Langmuir  circulation  increased  the  flux  of  gases  down  into  the  water 
column.  Rain  is  thus  able  to  effect  directly  the  surface  roughness  and  stress,  and  to  ef- 
fect indirectly  the  flux  of  gas  into  water  column. 

B.     LITERATURE  REVIEW 

1.     Surface  Effects  Due  to  Rain 

Reynolds  (1900)  observed  water  drops  falling  into  a  glass  vessel  and  he  postu- 
lated the  first  explanation  for  wave  attenuation.  From  his  study  of  vortex  rings,  he  be- 
heved  that  the  '"  knocking  down  "  of  waves  resulted  from  mixing  of  subsurface  water  by 
vortex  rings  formed  by  the  raindrops  hitting  the  surface.  As  the  vortex  ring  descended, 
it  entrained  fluid  from  the  moving  surface  layer  and  carried  this  fluid  downward. 
Reynolds  reported  that  the  nonmoving  fluid  from  below  would  rise  up  and  replace  the 
surface  fluid,  diminishing  the  motion  of  the  surface  which  is  essential  to  the  continuity 
of  the  wave. 

Later.  Ashton  and  O'Sullivan  (19-49).  Barnaby  (19-49).  and  Sainsbur>-  and 
Cheeseman  (1950)  further  investigated  Reynolds'  theory.  Each  of  these  groups  con- 
ducted experiments  involving  drops  of  dye  striking  a  water  surface.  Their  conclusions 
v^re  inconclusive  in  testing  Reynolds'  theory.  Although  vortex  rings  were  observed  in 
all  experiments,  the  rings'  effect  on  wave  attenuation  was  not  clearly  estabUshed. 

The  interaction  of  water  drops  on  surface  waves  was  explored  subsequently  by 
Siscoe  and  Levin  (1971).  They  found  two  ver\'  different  processes,  depending  upon  the 
raindrop's  point  of  incidence  on  the  wave:  a  reflection  process  and  an  absorption  proc- 
ess.   In  an  absorption  event,  the  drop  strikes  at  or  near  the  wave  crest  and  TKL  is  in- 


corporated  into  the  wave  field.  A  reflection  event  occurs  in  the  event  of  rain  striking  the 
trough,  and  the  TKE  of  the  drop  goes  into  producing  subsurface  turbulent  mixing.  They 
found  that  during  an  absorption  event  the  energy  supplied  by  an  impacting  drop  ampH- 
fied  the  surface  wave  field  which  propagated  out  from  the  point  of  impact.  Drops  were 
colored  by  dye  to  investigate  further  wave  formation  and  subsurface  interaction.  The 
dye  did  not  penetrate  the  surface  for  an  absorption  event,  while  during  reflection  events, 
dye  penetrated  to  the  subsurface  waters  due  to  drop-induced  mixing.  The  observed 
motions  of  the  dye  that  accounted  for  this  mass  transfer  were  vortex  rings  and  weakly 
defined  circulation  cells  that  reached  a  depth  of  5-10  cm  with  an  occasional  ring  pene- 
trating deeper. 

2.     Characteristics  of  Vortex  Rings 

\\"hen  raindrops  strike  the  sea  surface  a  portion  of  their  energy  is  transferred 
into  the  water  column  in  the  form  of  vortex  rings.  Keedy  (1967)  conducted  an  extensive 
study  on  the  formation  of  vortex  rings  by  free  surface  interactions.  Me  allowed  water 
drops  of  various  radii  to  fall  from  heights  of  2.5  inches  or  less,  and  he  observed  ring 
formation  and  penetration.  Keedy  found  that  vortex  rings  penetrate  to  a  maximum 
depth  of  approximately  10  cm.  Also  noted  were  density  differences  between  the  drop 
and  the  base  solution  which  were  important  in  determining  the  maximum  penetration 
depth  of  the  vortex  ring.  For  a  drop  with  a  density  lower  than  the  density  of  the  base 
solution,  the  vortex  ring  expanded  in  diameter  more  rapidly  than  if  both  drop  and  base 
solutions  were  of  the  same  density. 

Ma.xworthy  (1972)  observed  the  propagation  of  vortex  rings  prodi'-.ed  using  a 
membrane  apparatus  in  an  experimental  tank.  He  found  that  as  the  rings  traveled  they 
continuously  entrained  surrounding  fiuid  and  grew  in  volume.  At  the  same  time,  the 
rings  were  losing  vorticity  into  a  wake  and  decreasing  in  velocity.  Maxworthy  also  ob- 
served two  different  processes  when  two  similar  rings,  one  ahead  of  the  other,  came  into 


contact.  If  both  rings  had  the  same  velocity,  the  rearward  one  would  distort  and  squeeze 
into  the  forward  ring  through  the  rear  wake  area.  The  rearward  ring  was  then  wrapped 
around  the  inner  part  of  the  forward  ring  to  form  a  single  larger  vortex.  If  the  velocity 
of  the  rearward  ring  was  larger  than  that  of  the  forv^ard  ring,  the  rearward  ring's  entr\" 
into  the  forward  ring  caused  the  combined  ring  to  oscillate  longitudinally.  If  the  trailing 
ring's  velocity  is  much  greater  than  that  of  the  front  ring  Maxworthy  found  that  the 
newly  formed  ring  would  disintegrate  with  a  faster  ring  propagating  ahead  of  a  slower 


movmg  nng. 


Batchelor  ( 1967)  had  previously  found  a  similar  occurrence  involving  two  similar 
vortex  rings  and  suggested  an  interaction  analogous  to  "  leap  frog.  "  The  velocity  field 
of  the  rear  ring  has  a  radially  outward  component  at  the  position  of  the  forward  ring. 
As  the  radius  of  the  front  ring  increases,  travel  speed  decreases.  At  the  same  time,  the 
rearward  ring  undergoes  a  relative  increase  in  speed  and  passes  through  the  larger 
vortex,  becoming  the  lead  ring.  This  maneuver  is  then  repeated.  The  interaction  of 
these  rings  suggests  that  random  rings  produced  in  nature  will  collide,  changing  the 
rings  characteristics  and  velocity.  The  gross  eflect  will  be  to  reduce  the  penetration 
depth  expected  from  a  single  vortex  ring. 
3.     Rain-induced  MLxed  Layer  Changes 

Rather  than  study  individual  drops  and  vortex  rings,  Green  and  Houk  (1979) 
studied  populations  of  raindrops  and  the  subsequent  mixed  layer  formation.  Observa- 
tions of  changes  to  tbe  mixed  layer  depth  (MLD)  were  made  by  determining  the  location 
of  the  thermoclinc.  This  circumvented  the  need  to  consider  individual  vortices  in  favor 
of  an  integrated  turbulent  mixed  layer.  Green  and  Houk  found  that  a  mixed  layer  had 
developed  after  2(»  minutes  of  rain,  with  a  deeper  MLD  associated  with  the  larger  drop 
sizes.  The  depth  of  the  mixed  layer  in  salt  water  was  approximately  one  third  that  of 
freshwater  due  to  the  larger  downward  buoyancy  flux.    All  other  parameters  were  held 


constant.  Green  and  Houk  determined  that  the  majority  of  the  mechanical  energy  flux 
went  directly  into  subsurface  turbulence.  The  overall  effect  was  a  deepening  of  the  MLD 
as  the  energy  from  raindrops  diffused  down  into  the  water  column,  deepening  the  MLD 
10-30  cm  in  freshwater  and  5  to  10  cm  in  saltwater. 

4.     Manton's  Anal}1ical  Approach 
Manton  (1973)  presented  an  analytical  approach  to  explain  rain-induced  vortex  ring 
mixing  and  the  resulting  wave  attenuation.      Through  a   scale  analysis,  using  the 
raindrop's  terminal  velocity  (V)  and  the  horizontal  wind  speed  (U),  Manton  defined  a 
vertical  and  horizontal  eddv  viscositv: 


vv=/''iv  (2.i; 


Xn=f"Vl  (2.2) 

where /represents  the  volume  fraction  of  drops  in  the  air;  and  /represents  the  turbulent 
length  scale  deternuned  by  the  average  disiaiice  between  drops  in  the  air.  In  order  to 
use  vertical  and  horizontal  eddy  viscosities  in  the  momentum  equations.  Manton  needed 
to  prove  that  equations  2.1  and  2.2  were  independent  of  depth.  Through  analytical 
means  he  proved  that  the  shear  layer  thickness  for  a  wave,  on  the  order  (  -^  \\  was 
smaller  than  the  penetration  depth  of  vortices  (10-20  cm).  Erroneously,  the  shear  layer 
thickness  is  the  depth  associated  with  friction  from  a  boundary  and  not  the  depth  asso- 
ciated with  the  wave  energy,  1  k.  By  assuming  that  the  m-.ximum  mixing  depth  is  related 
to  vortex  ring  penetration,  mixing  should  only  occur  at  a  depth  of  10  to  20  cm.  which 
is  much  shallower  than  the  depth  associated  with  the  wave  energy.  Manton  solved  the 
linearized  momentum  equations  using  a  potential  velocity  solution,  implicitly  assuming 
irrotational  flow.  However,  when  mixing  is  present  the  flow  becomes  rotational.  It  is 
still  possible  to  predict  the  viscous  decay  of  the  waves  by  assuming  that  the  mixing  is 


uniform  over  depth  (Lighthill,  1978).  If  \\  and  v,,  are  assumed  to  be  equal,  Manton's 
results  would  agree  with  such  theory;  however,  Manton  used  separate  values  of 
v^  and  Vy,  to  incorporate  the  influence  of  wind  into  his  solution.  This  is  an  inappropriate 
use  of  the  viscous  decay  terms  in  the  momentum  equations.  For  a  vortex  ring,  the 
horizontal  and  vertical  mixing  scales  are  of  the  same  order,  thus  \\.  and  v^,  should  be  re- 
placed by  a  single  eddy  viscosity  coefficient,  v  ,  for  the  mixed  layer.  Wind,  which  is  not 
considered  in  this  study,  is  normally  introduced  as  an  additional  surface  stress  source 
term  in  the  momentum  equations  and  does  not  appear  as  part  of  the  viscous  dissipation 

term.   Thus  there  are  two  serious  objections  to  Manton's  effort: 

1.  introducing  mixing  over  the  entire  water  column,  and 

2.  incorporating  the  influence  of  wind  into  v. 

This  study  uses  a  numerical  solution  of  the  Navier-Stokes  equations  for  a  surface  gravity 
wave.  A  single  viscous  coefficient  v  will  be  applied  to  a  tliin  "  mixed  layer  "  and  the 
resulting  attenuation  will  be  estimated  and  compared  to  classical  theory  (Lighthill.  197S). 


III.     METHODS 

A.     THE  MODEL 

A  numerical  model,  SOLA-VOF  (Nichols  et  al,  1980),  developed  at  the  Los  Alamos 
National  Laboraton'  in  New  Mexico  was  used  for  this  study.  This  model  is  a  general 
two-dimensional  finite  differencing  Eulerian  code  which  allows  for  complicated  compu- 
tational domains  with  variable  boundar}'  conditions,  free  surface  shapes,  variable  surface 
tension  and  compressibility.  Behavior  of  the  fluid  is  determined  by  solving  the  Navier- 
Stokes  equations; 
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where  u  is  the  horizontal  velocity,  v  is  the  vertical  velocity,  g  is  the  body  acceleration  due 
to  gravity,  v  is  the  coefikient  of  kinematic  viscosity,  and  p  is  the  fluid  density.  In  addi- 
tion to  the  momentum  equations  the  incompressible  mass  conservation  equation: 
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ex 
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must  be  satisfied  at  each  time  step.  1  o  advance  a  solution  through  a  one  time  increment 

the  model  must  go  through  the  following  three  steps: 

\.  Explicit  approximations  of  Eq.  (3-1)  and  Eq.  (3-2)  are  used  to  compute  the  first 
guess  for  the  new  time-level  velocities  using  the  initial  conditions  or  the  previous 
time-level  values  for  all  advective,  pressure,  and  viscous  accelerations. 

2.  To  sati<;fy  the  continuity  equation.  Eq.  (3-3),  pressures  arc  itcratively  adjusted  in 
each  cell  and  the  \cIocity  changes  induced  by  each  pressure  change  arc  added  to 
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the  velocities  computed  in  step  (1).  A  repetitive  process  is  required  because  the 
change  in  pressure  needed  in  one  cell  to  satisfy  Ilq.  (3-3)  will  upset  the  balance  in 
the  four  adjacent  cells. 

3.  The  cell  fluid  function  defining  fluid  regions  is  then  updated  to  give  the  new  fluid 
configuration. 

In  this  manner  the  fluid  solution  can  be  advanced  through  all  desired  time  intervals. 
Horizontal  and  vertical  velocity  (w,v),  pressure  (P),  and  fractional  volume  of  fluid  in  a 
cell  (F)  are  calculated  at  each  time  step  at  the  cell  positions  shown  in  Fig.  2. 


-' 

\\*i 

^ 

. 

"i-i.\    ■ 
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■  "Ki.i 
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d                           ftv 

*                ox  1    ' 

Figure  2.  Typical  cell  configuration:  U  and  V  velocities  are  found  on  the  sides 
and  top  of  the  cell  with  pressure  and  fluid  percentage  of  each  cell  found 
in  tiie  center  of  the  cell. 
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The  fractional  volume  of  fluid  function  (F)  is  used  to  identify  grid  cells  that  contain 
fluid,  and  to  delineate  surface,  interior  or  boundar>'  cells.  A  free  surface  cell  is  defined 
as  a  cell  containing  a  nonzero  value  for  F  and  having  at  least  one  neighboring  cell  that 
contain  a  zero  value  for  F.  Cells  with  nonzero  F  values  and  no  empty  neighbors  are 
treated  as  cells  full  of  fluid. 

Introduction  of  variable  viscosity  into  the  code  required  a  change  to  the  code.  This 
was  accomphshed  by:  (1)  estabhshing  a  viscosity  grid  and  (2)  encoding  this  viscosity  grid 
into  the  viscous  terms.  The  viscosity  grid  was  established  by  identifying  three  areas,  v,, 
the  mixing  layer  at  the  surface  of  the  wave,  Vj,  a  transition  zone  from  the  mixed  layer  to 

V,  +  Vj 

the  turbulent  free  layer  defined  as  — - —  and  Vj,  the  base  layer,  (effected  by  molecular 
viscosity  only).  The  turbulent  mixing  layer  was  defined  as  those  cells  which  were  be- 
tween the  free  surface  and  the  desired  mi.xing  depth.  An  example  is  given  in  Fig.  ."^  de- 
picting the  mixing  layer  following  the  surface.  The  height  change  is  from  the  numerical 
initialization  of  the  wave. 


]?. 


Figure  3.  Kinematic  viscosity  grid  configuration:  The  setup  of  tlie  viscosity  grid 
where  v,  is  the  surface  mixing  layer,  Vj,  the  transition  layer,  and  Vj  the 
base  fluid  laver. 


The  viscosity  grid  was  redefined  at  each  time  step  resulting  in  a  moving  mixing  layer 
having  a  constant  thickness.  A  four  corner  average  was  then  used  to  obtain  a  value  of 
V  for  each  computational  cell.  The  single  derivative  form  of  the  momentum  equations 
(Harlow  1980): 
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was  used  to  reduce  numerical  noise  associated  with  calculating  second  derivatives  at  the 
boundan.'  of  the  two  layers. 

B.     THE  WAVE 

This  study  focused  on  the  attenuation  of  deep  water  waves  resulting  from  rain- 
induced  mixing.  To  simulate  this  a  computational  domain  was  initialized  representing 
a  fluid  vessel  or  tank  1.4  m  wide  and  3.0  m  deep.  A  Mean  Fluid  Height  (MFH)  of  2.8 
m  was  selected.  The  tank  depth  of  one  wavelength  was  chosen  to  ensure  that  bottom 
effects  would  be  negligible  during  the  propagation  of  the  wave.  According  to  deep  water 
wave  theory,  velocity  at  the  base  of  the  tank  is  smaller  than  the  surface  velocities  by  a 
factor  of  ^-'''"  =  (?--''<  1  .  Furthermore,  a  small  amplitude  wave  (3  and  6  cn;i)  was  chosen 
to  maintain  stability  of  the  wave.  During  the  early  part  of  the  project  the  wave  was 
found  to  break  if  the  amplitude  became  greaier  than  10  cm  due  to  the  manner  in  which 
the  code  assigned  velocities  at  the  free  surface.  Normally  a  surface  velocity  is  set  to  be 
the  same  value  as  the  nearest  horizontal  neiglibor  cell.  Depending  on  the  direction  of 
wave  propagation  this  horizontal  movement  of  velocities  occasionally  produced  an  un- 
reahstic  velocity  discontinuity  in  the  vertical.  This  problem  was  corrected  by  making  the 
nearest  neighbor  cell  directly  below  the  surface  cell. 

A  2.S  m  wavelength  standing  wave  was  chosen  to  minimize  continuation  boundar>' 
problems  associated  with  a  progressive  wave,  and  to  facilitate  a  simple  initial  condition 
having  no  velocities.  A  standing  wave  is  actually  two  progressive  waves  propagating  in 
opposite  directions  given  by:  ?/  =  ;/„  cos  kx  cos  w/  where  7/^,  is  the  wave  amplitude,  k  is 
the  wave  number,  and  o  is  the  radial  frequency.    A  progressive  wave  was  simulated  by 


initializing  one-half  wavelength  of  the  standing  wave  using  the  instantaneous  wave 
height:  >/  =  MFH  —  ij^  cos  kx  at  t  =  0  for  a  domain  ofO<A-<— .  The  waves  initial 
jnfiguration  is  seen  in  Fig.  4  under  steady  state  conditions.  Velocity  vectors  are  also 
plotted  for  the  initial  steady  state  condition  and  are  represented  by  a  vector  of  zero 
len2th. 
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Figure  4.      Initial  \\a\e  coiiliguration  for  all  experiments:     1  he  one-lialf  wavelength 
at  rest  prior  to  being  released. 

The  wave  was  released  from  rest  and  observed  as  it  propagated  through  the  lank.  \\'a\e 
dynamics  were  reviewed  and  compared  to  theoretical  values  to  ensure  that  the  wave  and 
code  had  been  properly  developed. 

Validation  of  the  wave  was  accomplished  by  examining  the  period,  the  vertical  U 
velocity  profile,  the  potential  and  kinetic  energy,  and  the  attenuation  characteristics  of 
the  wave.  According  to  deep  water  wave  theory  a  period  of  1.34  s  is  derived,  in 
Fig.  5  and  Fig.  6  the  wave  is  seen  to  take  approximately  1.3  s  to  propagate  from  a 
maximum  amplitude  at  the  right  wall  to  the  left  and  back  to  the  right  wall  coinciding 
with  theorv. 
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Figure  5.      Wave  propagalion  through  the  tank:     The  wave  is  observed  propagating 
from  rest  toward  the  left  wall  of  the  tank. 


16 


T  =  0.90  a   CYCLE  =  90 


T=  1.00  9   CYCLE  =  100 


T  =  1.10  s   CYCLE  =  110 


T  =  1.20  s   CYCLE  =  120 


T  =  1  30  s   CYCLE  =  130 


.01  m/s 


T  =  1.40  3   CYCLE  =  140 


rt 


Figure  6.      Wave  propagation:     The  wave  is  now  propagating  toward  the  rigiit  wall 
of  the  tank. 
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Figure  7.      Natural  log  of  u  vs  depth. 

To  validate  the  velocity  field  a  vertical  U  velocity  profile  was  examined  and  compared 
to  theor\':  u  =  ii^  e-'".  The  theoretical  velocity  decay  profile  was  derived  and  plotted, 
Fig.  7.  The  numerical  profile  was  consistent  with  theory  except  for  an  inconsistency 
seen  near  the  bottom  of  the  tank  which  may  be  due  to  numerical  noise  contamination. 

The  energy  of  a  standing  wave  is  known  to  be  periodic  alternating  from  Potential 
Energy  (PE)  to  Kinetic  Energy  (KE)  over  time.  Theoretical  values  of  the  total  energy 
density  are  obtained  by  adding  the  energy  density  of  the  two  progressive  waves, 
( '/2  p  g  ijI  ^.)  each  with  half  the  amplitude,  (;/^  =  Yi  >/o)  of  the  standing  wave.  For  a  wave 
amplitude  of,  ijq  ,  6  cm  the  total  energy  is  V^pg^il  =  12.348  J  and  for  a  3  cm  amplitude, 
3.US7  J.   To  compute  energy  numerically,  PE  density  was  calculated  from: 
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PE=V2pgri^dxdz  (3.6) 

by  summing  over  the  surface  of  the  wave,  where  r\  is  found  to  be  the  distance  the  wave 
surface  is  above  or  below  the  MFH.   KE  density  was  found  by  summing: 

KE=  VipW^dxdydz  (3.7) 

over  the  domain  of  the  tank  where  V^  =  (i^  +  v^)  for  each  cell.  The  energy  for  a  6  cm 
amplitude  wave  propagating  under  normal  conditions  of  molecular  viscosity 
(v  =  10"''  w?  5-')  is  seen  in  Fig.  8.  The  peak  energy  values  are  consistent  with  theory. 
Initially  PE  is  a  maximum  and  KE  is  zero.  As  the  wave  propagates  PE  approaches  zero 
and  KE  reaches  a  maximum  at  each  one-quarter  cycle.  PE  becomes  a  maximum  at  each 
one-half  cycle  as  KE  approaches  zero. 
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Figure  8.      Potential  and  Kinetic  energy  vs  time:     7  lie  sinusoidal  nature  of  PE  and 
KE  is  depicted  for  the  wave  as  it  propagates  over  time,  v  =   10"*/?i^5"'. 

Kinematic  viscosity  was  used  to  represent  the  turbulent  mixing  occurring  in  the 
water  column.  The  values  used  ranged  from  10'  ni^  5-'  representing  molecular  effects  to 
10-^/7?' 5-'  representhig  high  turbulent  conditions.  This  range  included  values  found  by 
Manton:  3  x  lO-^m^s"'  for  a  light  rain  event  and  1.6  x  10'^  m^  s~^  for  a  hea\7  rain 
event.  This  range  also  included  Green  and  Houk's  estimate  for  the  mixed  layer  eddy 
viscosity  coefTicient,  5  x  10-'  m^  5-'.  The  kinematic  viscosity  was  kept  constant  over  the 
entire  water  column  to  determine  if  attenuation  effects  were  following  theory.  Both 
amplitudes  were  used  for  this  test.  1  wo  different  methods  were  used  to  determine  the 
rate  of  energy  change  over  time.   The  slope  of  the  decay  of  peak  kinetic  energy  was  used 
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as  a  direct  representation  of  the  rate  of  energy  change.  From  Fig.  9  the  attenuation  due 
to  V  =  10'^  ni?  5-'  gives  a  normalized  slope  of  0.15  which  is  nearly  consistent  with  theor>", 
Avk^  =  0.2.  A  more  accurate  method  was  found  using  the  internal  dissipation  function 
(Schlichting.  1960): 
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This  represents  the  decrease  of  energy  from  viscous  effects  and  does  not  include  the  re- 
distribution of  energy  into  the  wave.  At  each  time  step  the  change  of  energy  was  cal- 
culated and  normalized  by  the  instantaneous  kinetic  energy.  Final  processing  of  the  data 
was  then  conducted  by  taking  the  average  of  the  rate  of  change  of  energy  divided  by 
KE,  i.e.,<l)/A'£.  for  KE  >  1.  This  was  required  due  to  an  erroneous  spike  developing  as 
KE  became  less  than  one.  as  seen  in  Fig.  10.  This  was  done  for  both  the  3  and  6  cm 
amplitude  cases,  the  results  of  which  are  gi\en  in  Tables  1  and  2.  showing  good  agree- 
ment with  theor>-. 
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Figure  9.  Potential  and  Kinetic  energy  vs  time:  The  sinusoidal  nature  of  PE  and 
KE  are  depicted  for  the  wave  as  it  propagates  over  time, 
V  =   10-2A7J-5-'. 

C.  THE  MIXED  LAYER  PROBLEM 

Due  to  observations  by  Keedy,  Siscoe  and  Levin,  and  Houk  and  Green  two  mixing 
depths,  10  ano  20  cm  were  used  in  this  study.  Ihe  10  cm  depth  was  r  presentative  of 
the  maximum  penetration  depth  of  vortex  rings  formed  by  rahidrops  striking  the  ocean 
surface  while  20  cm  represented  an  increased  mixing  depth  penetration,  but  occurring 
more  infrequently.  Kinematic  viscosity  of  the  base  layer  was  set  at  10"*  m^  5"'.  'Ihe 
surface  mixing  layer  was  varied  from  10~^  to  10^^  ah' 5-'.  As  for  the  single  layer  case, 
<\>IKE  was  calculated  at  each  time  step  and  averaged  for  all  values  of  KE  greater  than 
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Figure   10.       Ra(e  of  energy  change:     Attenuation  for  cl)  (— )  and  <]^IKE  (---)  with 
V  =   10-2  m^  5-'. 

one.  These  results  were  compared  to  the  expected  theoretical  values  of  the  two  layers. 
The  theoretical  values  shown  are  -^  where  v,  is  the  mixing  if  the  viscosity  coeflicient  in 
the  mixed  layer  is  used  and  Vj  is  the  mixing  if  the  viscosity  coefficient  (molecular)  of  the 
lower  layer  is  used. 

The  effects  of  increasing  v  in  the  top  10  cm  layer  is  expected  to  increase  attenuation, 
although  less  than  expected  if  the  entire  water  column  is  more  turbulent.  Doubling  the 
mixing  layer  to  20  cm  is  hypothesized  to  double  the  attenuation  seen  in  the  10  cm  sur- 
face layer  case.  By  increasing  the  penetration  depth  more  kinetic  energy  should  be  taken 
out  of  the  wave,  i.e., dampen  the  wave. 
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The  results  shown  in  Tables  3  and  4  are  for  a  10  cm  mixed  layer  and  in  Table  5  for 
a  20  cm  mixed  laver. 


Table  3.     RESULTS   FOR   A   3   CM    AMPLITUDE  WAVE,    10  CM   MIXING 
DEPTH. 


Xu 

(D/A'£ 

Theory  Svk-  vjv^ 

l.E-4 

.00162 

.0040  .00004 

3.E-4 

.00477 

.0120.00004 

l.E-3 

.0157 

.0400  .00004 

1.6E-3 

.0243 

.0640,  .00004 

l.E-2 

.1269 

.4000  .00004 

Table  4.     RESULTS   FOR  A   6  CM   AMPLITUDE   WAVE,    10   CM   MIXING 
DEPTH. 


Nu 

CD/AT 

Theor\'  8\7c*  \\l\\ 

l.E-4 

.00 1 65 

.0040  .00004 

3.E-4 

.00489 

.0120  .00004 

l.E-3 

.0159 

.0400,  .00()04 

1.6E-3 

.0252 

.0640  .00004 

1.E-: 

.1379 

.4000  .000(:)4 

Table  5.     RESULTS    FOR   A   6   CM    AMPLITUDE   WAVE,   20   CM    MIXING 
DEPTH. 


Nu 

cD/A"£ 

Theory-  SvA-  v,/'v-, 

l.E-4 

.00259 

.0040  .00004 

3.E-4 

.00773 

.0120  .00004 

l.E-3 

.0256 

,0400  .00004 

1.6E-3 

.0409 

.0640  .00004 

l.E-2 

.2440 

.4000  .00004 

1.     Single  Layer  Wave 

Using  molecular  mixing,  v  =  10"*  m^  5-'  ,  (Tables  1  and  2)  a  3  cm  amplitude  wave 
required  46,000  wave  periods  and  a  6  cm  amplitude  wave  required  82,000  wave  periods 
to  reduce  KE  to  an  arbitrary"  cutoff  value  of  0.5  N-m.  By  comparison,  the  rate  of  energy 
change  for  a  3  cm  amplitude  wave  with  a  10  cm  mi.xed  layer  {\\  =  10-  ni-  s'^)  showed  a 
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3,000  fold  increase  in  the  rate  of  energy  loss,  requiring  only  14  wave  periods  to  dampen 
the  wave.  An  almost  identical  rate  of  energy  loss  is  seen  for  the  6  cm  amplitude  wave 
with  dampening  occurring  in  22  wave  periods.  If  the  mixing  depth  was  doubled  to  20 
cm  there  was  a  6,200  fold  increase  in  the  rate  of  energy  change  which  is  nearly  double 
that  of  the  6  cm  amplitude  wave  with  a  10  cm  mixing  depth.  Thus  by  increasing  tur- 
bulent mixing  in  a  thin  surface  layer  an  increase  in  the  rate  of  energy  change  (3000-6000 
fold  or  greater)  can  be  produced.  This  is  especially  noteworthy  for  short  gravity  and 
capillar}-  waves  whose  energy  is  mostly  contained  within  the  mixing  layer. 
2.     Rain's  Effects  on  Waves 

a.  10  cm  Mixing  Depth 

To  model  the  effects  of  rain,  three  viscosities  were  directly  associated  with 
rain  events: 

1.  3.0  X  lO"*  m^  5~'  for  a  hght  rain  event   (5  mmlhr  rain  rate)  (Manton) 

2.  1.6  X  1()-^  m-  5"'  for  a  hca\7  rain  event   (12  mmlhr  rain  rate)  (Manton) 

3.  5.0  X  10"-*  m-  5"'  for  a  hea\y  rain  event   (14  mmlhr  rain  rate)  (Green  and  Ilouk) 

Viscosities  (1)  and  (2)  were  used  in  th'S  study.  Comparing  the  results  for  these  two 
viscosities  showed  that  as  the  rain  rate  and  the  turbulent  mixing  increased,  an  increase 
in  the  attenuation  of  the  wave  also  occurred.  The  time  interval  required  to  dampen  a  3 
cm  amplitude  wave  with  lU  cm  mixing  depth  was  380  wave  periods  for  light  rain  and  80 
wave  periods  for  hcaNy  rain.  Doubling  this  wave  amplitude  and  maintaining  a  constant 
10  cm  mixing  depth,  is  found  to  increase  the  wave  periods  required  to  dampen  the  wave. 
During  a  light  rain  event  the  wave  is  dampened  within  655  wave  periods  and  126  wave 
periods  during  a  hea\y  rain. 

b.  20  cm  Mixing  Depth 

As  shown  in  Table  5,  doubling  the  mixing  depth  did  not  double  the  rate  of 
energy  change  as  hypothesized,  although  a  significant  (50-80°'o)  change  did  occur.  For 
a  wave  of  2.8  m  wavelength  and  6  cm  amphtude,  a  hght  rain  event  caused  the  wave  to 
be  dampened  within  420  wave  periods  and  .dthin  80  w  ve  periods  for  a  hea\7  rain  event. 
This  is  a  reduction  of  approximately  38*^/'o  of  the  wave  periods  required  to  dampen  the 
same  wave  with  a  10  cm  mixing  depth. 
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IV.     CONCLUSION 

Strong  rain-induced  mixing  within  a  ver>'  thin  surface  layer  will  increase  the  atten- 
uation rate  of  surface  gravity  waves  by  several  orders  of  magnitude  over  the  rate  of  at- 
tenuation due  to  molecular  viscosity  alone.  In  this  case  study  a  single  wavelength  (2.8m) 
was  used  together  with  two  mixed  layer  depths  (10  and  20  cm).  In  the  case  of  strong 
mixing  (v  =  10'^  m^  s'^)  with  a  10  cm  thick  mixed  layer,  the  wave  was  dampened  with 
12-24  wave  periods  compared  to  46,000  wave  periods  for  molecular  viscosity  alone. 
Using  estimates  of  the  eddy  mixing  coefficients  for  hea\7  and  light  rain  from  Manton 
(v  =  \.6x\0-^  m^  s'^  and  v  =  3.TlO"^mV)  .  the  wave  is  still  dampened  within  126  and 
655  wave  periods,  respectively.  Increasing  tlie  mixed  layer  depth  to  20  cm  produced  a 
significant  increase  in  attenuation,  reducing  the  nunber  of  wave  periods  required  to 
dampen  the  wave  to  80  and  420  wave  periods,  respectively  for  heav\'  and  and  Ught  rain. 
This  influence  of  rain  on  wave  attenuation  has  several  potential  applications. 

Scatterometer  and  SAR  observations  of  the  sea  surface  depend  on  the  presence  of 
short  gravity  and  capillar}'  waves.  These  waves  are  generally  shorter  than  the  case  study 
wavelength  (2.S  m)  and  they  will  be  even  more  strongly  attenuated  by  rain-inducid 
mixed  layers.  Thus,  this  theor>"  can  be  used  to  explain  the  observed  SAR  imagery  in 
rainy  regions  (dark  image)  and  may  allow  a  correction  to  be  applied  to  scatterometer 
algorithms  used  to  measure  wind  speed  from  small  scale  ocean  roughness.  This  should 
increase  the  accuracy  of  weather  forecasts  lor  marine  shipping.  The  damping  of  slightly 
longer  gravity  waves  present  in  ship  wakes  is  also  of  interest.  Rain  will  attenuate  these 
waves  and  thus  mask  the  ship's  wake,  decreasing  the  chance  of  detection  by  other  re- 
mote sensing  techiques.  Rainfall  measurements  can  also  be  incorporated  into  wave  and 
surf  forecasts  as  wave  length  dependent  attenuation.  This  should  increase  the  accuracy 
of  such  models  by  decreasing  the  >;   ergy  contained  l.  the  high  frequency  spectrum. 

This  study  has  developed  a  boundar\'  layer  that  resembles  the  shape  of  the  water 
surface  which  may  assist  in  simulating  mixing  problems  more  reahstically.  The  eflect 
of  turbulent  kinetic  energy  could  be  introduced  numerically  by  using  the  turbulent 
kinetic  energy  equation.  This  equation  would  provide  a  "Reynolds'  eddy  viscosity  coef- 
ficient" for  each  cell  which  could  then  be  used  directly  rather  than  by  defining  an  arbi- 
trarv  mixed  lavcr.   This  will  allow  variable  surface  cravit\  wave  attenuation  as  the  mixed 
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layer  deepens.   Through  this  method  the  efTect  of  a  finite  duration  rain  squall  on  mixed 
layer  depth  and  wave  attenuation  might  be  estimated. 

This  study  focused  on  one  wavelength,  two  amplitudes  and  two  mixing  depths.  It 
is  recommended  that  further  numerical  studies  be  conducted  with  larger  wavelengths, 
higher  wave  ampHtudes  and  deeper  mixing  depths.  These  factors  will  identify  cutoff  re- 
gions where  rain-induced  attenuation  no  longer  has  an  effect  on  the  wave.  Wave  tank 
experiments  would  also  benefit  numerical  modeling  to  assist  in  parameter  formations  for 
the  many  inputs  necessar>"  for  turbulent  kinetic  energy  models. 
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